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Influence of tetragonal distortion on the magnetic and electronic properties 
of the Heusler compound Co2~TiSn from first principles 
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Using the full potential linearized augmented plane wave plus local orbitals method we determine ab-initio the 
lattice parameters of tetragonally distorted Co2TiSn in the L2i structure. The tetragonal lattice parameter 
c is determined as a function of the lattice parameter a by energy minimization. The change in total energy 
is found to be only a few fc#T with respect to room temperature. The spin polarizations as well as the 
magnetizations are stable against small lattice distortions. It is shown, that the volume is not constant upon 
distortion and that the volume change is related with significant changes in the magnetization and the gap 
energy. 

PACS numbers: 61.50.Ah, 71.20.Be, 75.50.Cc 
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C02 based Heusler compounds have been subject of 
extensive studies in the context of spintronics during the 
last decade. They are particularly of interest because 
of their predicted full spin polarization at the Fermi 
edge, making them so-called ferromagnetic half-metals. 
Fully epitaxial magnetic tunnel junctions grown on Cr 
buffered MgO(OOl) with a C02 based Heusler compound 
as one electrode can exhibit more than 200 % tunnel- 
magnetoresistance ratio (TMR) at room temperature 
(RT); e.g., 386% with C02FeAlo.5Sio.3i or 216% with 
Co2MnSi 2 were demonstrated. Fully epitaxial junctions 
without Cr buffer exhibiting high TMR of 170% at RT 
and two Co2MnSi electrodes,— as well as junctions with 
an epitaxial Co 2 MnGe electrode with 160 % TMR at RT 
are reported^ The {011} lattice spacing of MgO differs 
by a few percent from the typical Heusler {001} spacings. 
Thus, one expects to find significant strain in thin epi- 
taxial films of these compounds on MgO if they are not 
buffered by a material with intermediate lattice parame- 
ter (as Cr). Due to symmetry, a tetragonal distortion is 
the most likely one. By assuming a constant volume of 
the unit cell, Block et al. have calculated the spin polar- 
ization of Co2CrAl under tetragonal distortions.— They 
find that the spin polarization will remain stable if the 
distortions are of the order of 2-3%. 

Here, we present band structure calculations of 
Co2TiSn (CTS) in the L2i structure under tetragonal 
distortions without the assumption of a constant unit 
cell volume. There are numerous theoretical studies of 
CTS using diverse band structure schemes.— ~— The re- 
sults with respect to magnetization and half-metallicity 
differ considerably with the applied method. Kandpal 
et al. and Miura et al. predict CTS to be a promising 
low-moment (2 /x#) ferromagnetic half-metal with a high 
stability against the Co-Ti swap disorder . 11 i 12 

For this study the Elk package,— an all-electron full- 
potential linearised augmented-plane wave plus local or- 
bitals (FP-LAPW+lo) code, was used. All calculations 



were performed with a k-point grid of 16 x 16 x 16 fc-points 
on the irreducible wedge of the Brillouin zone, yielding 
349 fc-points. In all calculations the primitive tetrago- 
nal cell was used. A smearing width of 0.027 eV was set 
and the number of initial empty states was fixed to 40. 
For the exchange-correlation the generalized gradient ap- 
proximation (GGA) was employed in the Perdcw-Burkc- 
Ernzerhof (PBE) implementation. 14 Spin-orbit coupling 
of the valence states was neglected, the muffin-tin radii 
were automatically optimized to be nearly touching. The 
total energy was converged to 2.7T0 -3 eV. The maximum 
length of the G + k vectors was given by the muffin-tin 
radius -R^in y i & ^min x max |G + k| = 7. For some criti- 
cal geometries the total energy and the magnetic moment 
were checked for convergence with respect to the cut-offs 
for the basis set; the parameters used for our calculations 
were found to provide well converged results. 

First, the cubic equilibrium lattice parameter a eq was 
determined by total energy minimization (Figure [1} . We 
find 6.12 A from theory, whereas the experimentally re- 
ported lattice parameter for bulk Co2TiSn is 6.0718A. 12 
This corresponds to an error of 0.8%. However, it is a 
common finding, that lattice parameters predicted from 
density functional theory energy minimizations using the 
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FIG. 1. Relaxation of the cubic unit cell with total energy 
(circles) and magnetic moment (squares). The range of lattice 
parameters corresponds to ± 10 % strain in 2 % steps. For the 
minimum determination more points were calculated, which 
are not shown here. 
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FIG. 2. Total energy (top) and magnetic moment (bottom) as 
a function of the tetragonal lattice parameter c, displayed for 
the different distortion parameters A. The energy is set to zero 
at the global energy minimum, i.e., at the total energy of the 
relaxed cubic configuration with a cq = 6.12 A. The magnetic 
moments corresponding to -E m i n (A,c) are tagged with large 
markers. 



GGA exchange-correlation tend to overestimate the real 
lattice parameter slightly, whereas local spin density ap- 
proximation (LSDA) calculations without gradient cor- 
rections tend to underestimate them?^ 

Biaxial strain was simulated by tetragonal lat- 
tice distortions using a = A ■ a eq with A = 
0.94, 0.96, . . . , 1.04, 1.06. The equilibrium value c(A) was 
then determined by total energy minimization with fixed 
a. For every A considered, the calculations were per- 
formed on 21 values of c, making a total of 147 data 
points. 

The total energies and the magnetic moments obtained 
from the calculations are displayed in Figure El It is 
immediately clear that the cubic (A = 1.00) configuration 
is the most stable state with minimum energy. 

The total energy for biaxial strain of ± 2% is increased 
by 25 - 30 meV, i.e., about the typical thermal energy at 
room temperature. Compared to the energies involved 
with homogenous lattice compression/expansion (Figure 
[1]), this energy scale is rather small. Further, the minima 
of the A ~ 1 configurations are relatively broad. This al- 
lows for further distortion of the lattice when additional 
influences get involved. In the case of a thin film growing 
on a single crystalline substrate, these are, among others, 



the surface bonding of the substrate and the relaxation 
of the film at its surface (and, in case of sputtering, in- 
teractions with the sputter gas). In fact, the structure 
has to be considered as very soft with respect to small 
non-relaxed tetragonal distortions, allowing thin films of 
CTS to attain very different real structures in dependence 
on their growth conditions, post-growth treatments and 
long-range structure. According to our results, a change 
in the lattice parameters of a CTS thin film with vary- 
ing growth conditions or annealing temperature cannot 
solely be ascribed to chemical ordering and recrystalliza- 
tion, but also to strain caused by epitaxial matching with 
the substrate. 

The magnetic moments in dependence of c(A) reveal 
a complex response of the density of states and the ex- 
change interaction to the distortions. Here, we take a 
magnetic moment of 2 fis / f-u. as indicative for a half- 
metallic state, which is confirmed by the density of states 
(Figure [3]). The most striking aspect here is the dif- 
ferent behaviour if the lattice is epitaxially compressed 
(A < 1) or expanded (A > 1). With increasing compres- 
sion the half-mctallicity vanishes quickly or even mag- 
netism is lost completely (A = 0.94). In the case of epi- 
taxial expansion the situation is different, half-metallicity 
is conserved here. The discontinuous transition from the 
ferromagnetic to the non-magnetic case with decreasing c 
pins the system in the ferromagnetic, half-metallic state. 
The reason for the loss of the full magnetic moment for 
the most non-relaxed geometries considered here is a clos- 
ing of the gap, which is due to the breaking of the cubic 
crystal field symmetry. This in turn is responsible for the 
formation of common bonding and anti-bonding hybrid 
bands, which drive the gap formation. 8 The closing of 
the gap can not be ascribed to any atomic site in partic- 
ular; Co, Ti as well as Sn exhibit a non-zero DOS at the 
position of the gap when the magnetization is less than 
2/x B /f.u. 

A sharp transition from magnetic to non-magnetic is 
observed if c is reduced below a critical value c cr it , which 
depends weakly on A: in the compressed case, Cent is a 
little smaller than for A = 1, whereas it is a little bit 
larger in the expanded case. This is clear, considering 
that the geometry is closer to the cubic for a smaller c 
with A < 1 and for a larger c with A > 1 . It is remarkable 
that c cr it does not depend on the actual value of A (on the 
scale of our c(A) sampling). However, this critical lattice 
parameter of 5.5 - 5.8 A is also found for the homoge- 
nous compression of the lattice (Figure [l}, indicating a 
nearly geometry-independent minimum distance for the 
ferromagnetic coupling of the Co atoms. 

The comparison of distortion calculations with con- 
stant volume and the relaxed calculations (Figured]) re- 
veals that the assumption of a constant volume is only 
partially justified. For A < 1 it turns out, that the vol- 
ume remains essentially constant. In contrast, for A > 1 
the volume increases significantly. For A = 1.06 we ob- 
serve a volume increase of 6.8%. This large discrepancy 
is caused by the above mentioned energetic pinning to the 
ferromagnetic state. Here we also see the largest differ- 
ence between the constant-volume calculations and the 
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FIG. 3. Densities of states for all A-values with relaxed c. Spin-up is always presented on the left, spin-down on the right side 
of the plots. The Fermi energy is at zero energy. 



fully relaxed ones: whereas the constant-volume calcula- 
tions predict a very small window of A in which the com- 
pound is magnetic, the fully relaxed calculations rather 
predict a quite large stability of the ferromagnetism and 
even the half-metallic gap. 

In summary we have shown by ab-initio band structure 
calculations that the cubic crystal structure of Co2TiSn 
allows easily for tetragonal distortion, which we therefore 
expect to occur in epitaxially grown thin films. The half- 
metallicity is retained over a large range of distortion. 
We provide an additional explanation for the typically 
observed large spectrum of lattice parameters of Heusler 
compound thin films in dependence on deposition and 
post-treatment conditions. In comparison with calcula- 
tions that assume a constant volume, we show that this 
assumption is not justified due to the energetically favor- 
able ferromagnetism. Thus, when dealing with distor- 
tions in Heusler compounds, it is necessary to consider 
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relaxed geometries. We assume, that these findings are 
representative for other C02 based Heusler compounds, 
although the results may vary in detail. Further investi- 
gation on this subject is necessary. 
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FIG. 4. Top: Relaxed unit cell volume as a function of the 
distortion parameter A. Bottom: Total energy (circles) and 
magnetic moment (squares) for constant volume (open sym- 
bols) and relaxed (filled symbols) calculations. 



